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We study gravitational leasing by a recently proposed black hole solution in Loop Quantum 
Gravity. We highlight the fact that the quantum gravity corrections to the Schwarzschild metric 
in this model evade the ‘mass suppression’ effects (that the usual quantum gravity corrections are 
susceptible to) by virtue of one of the parameters in the model being dimensionless, which is unlike 
any other quantum gravity motivated parameter. Gravitational lensing in the strong and weak 
deflection regimes is studied and a sample consistency relation is presented which could serve as a test 
of this model. We discuss that though the consistency relation for this model is qualitatively similar 
to what would have been in Brans-Dicke, in general it can be a good discriminator between many 
alternative theories. Although the observational prospects do not seem to be very optimistic even 
for a galactic supermassive black hole case, time delay between relativistic images for billion solar 
mass black holes in other galaxies might be within reach of future relativistic lensing observations. 
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I. INTRODUCTION 

The construction of a consistent theory of quantum gravity describing our universe remains a cherished dream of 
theoretical physicists. Two foremost candidates for a theory of quantum gravity at present are String Theory and Loop 
Quantum Gravity (LQG) among many. There are multiple pointers to believe that study of black holes is likely to yield 
a few insights for revealing the quantum nature of gravity and hence the quantum nature of spacetime. The singularity 
problem, the entropy puzzle and the information loss paradox (especially in light of some recent developments [1-3]) 
leads one to believe that black hole might be fundamentally quantum object in some sense. Motivated by these 
developments, it has been recently argued that precision interferometry of lensed pulsar images may shed light on the 
quantum gravitational processes near black holes [ ]. Large enough black holes would be described by a classical black 
hole metric which can encode the quantum corrections in a suitable way. For example, in the spherically symmetric 
case, the deviation from the Schwarzschild metric can be characterized by a parameter (apart from mass) involved in 
the corresponding quantum gravity theory such as the string tension in the string theory case or a non-commutative 
parameter in certain non-commutative quantum geometry scenarios or the Barbero-Immirzi parameter for LQG. Of 
course one expects the quantum gravity effects to be relevant when spacetime curvature reaches the Planckian regime 
which is reflected by the fact that string tension or non-commutativity parameter have dimension of [L]^ and hence 

their effect on observables is mass suppressed i.e. typically the observables are corrected by factors of 1 -I- 
where Iqg is a quantum gravity scale and m is an arbitrary positive power. Hence the correction is important for 
very low mass black holes, indeed microscopic ones if Iqg is near Planck scale £p^. Nonetheless it has been of some 
interest to study the effects of such corrections both quantitatively and qualitatively for string inspired scenarios [5] 
and non-commutative black holes [6]. These studies have concentrated on gravitational lensing as an observational 
probe, especially in the strong field regime. The strong field lensing coefficients are interesting in their own right too, 
owing to their relationship with quasinormal modes [7] and absorption cross sections from black holes [8]. In this 
paper we study gravitational lensing in a loop quantum gravity inspired scenario [9] both in the strong and weak field 
regimes. The peculiarity of the particular solution we study, as we shall discuss in detail later, is that the quantum 
gravity corrections, apart from entering into the metric function as a dimension term viz. the minimal area quanta. 
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also enter as a dimensionless parameter, viz., the ‘polymeric parameter’. This is in contrast to the stringy and non- 
commutative cases discussed before and makes the lensing observations evade the mass suppression effect and is by 
far the most observationally favorable quantum gravity scenario, though still quite difficult observational challenge in 
itself. 

The paper is organized as follows. In section Section II the loop black hole spacetime is introduced and in sections 
Section III and Section IV respectively, the weak field and strong field lensing signatures are studied. After a brief 
discussion of observational prospects in section Section V and the contrast with other quantum gravity scenarios in 
section Section VI, a sample consistency relation between the far field and relativistic Einstein Rings is presented in 
section Section VII. Finally we conclude with discussions in section Section VIII. 


II. A BLACK HOLE SOLUTION IN LQG 

Let us briefly discuss the black hole solution in LQG. LQG is a potential candidate theory of quantum gravity 
[10, 11]. In this theory, su{2) valued connection (evaluated for the holonomies along closed loops) and triad densities 
serve as the canonical variables for quantization. The solution of this theory is obtained in terms of constraint equations 
(Gauss, Diffeomorphism and Hamiltonian constraints respectively) whose solutions will describe the dynamics of space 
time. Though the theory is well understood at the kinematical level in terms of solutions of Gauss and diffeomorphism 
constraints known as ‘spin-networks’, solution for all the constraints in the full generality has yet not been achieved. 
The theory tentatively suggests discrete structure of space-time at the kinematical level. Geometric quantities like 
area, volume are shown to have discrete spectra (for the Euclidean case which can be carried over to Lorentzian case at 
the kinematical level). However, the full dynamical structure of the theory is missing as the solving the Hamiltonian 
constraint remains an open issue. 


A. Connection Formulation for Schwarzschild 


In Loop Quantum Gravity, black holes in different spacetimes have been studied [9] employing the semiclassical 
techniques in the minisuperspace quantization scheme. In this work we consider a black hole model for the quantum 
corrected Schwarzschild spacetime [12], which is characterized by a ‘polymeric parameter’ 5, which is a measure of 
‘quantum correction’ to the classical description. In order to study the loop corrected behavior of Schwarzschild black 
hole, one starts with the homogeneous but anisotropic Kantowski-Sachs spacetime expressed in terms of Ashtekar 
variables. 


A = CTsdx + bT2d9 — bri sin ddcj) + T 3 cos 9d(j), 

~ ^ 

E = PcT 3 Sintl— -|-Pbr 2 sm 6 '— 


( 1 ) 


with Ti € su{2). Therefore, the basic canonical variables in this approach are {b, c;pb,pc}- The Hamiltonian constraint 
in terms of these variables is given as 


Ch 


/ 


Ndx sin 9d9d(j) 
SttGnj'^ 


vW^\ 


+ 2bc , 


( 2 ) 


with N giving the lapse function on the spatial surface given by the the metric related to the triad Ef as 

= det(Q)q“^ 

and 7 being the “Barbero-Immirzi” parameter which appears in the action of the theory, which at the classical level, 
is an equivalent canonical formulation of Einstein-Hilbert action. This parameter remains free in the quantization 
scheme and characterizes different unitarily inequivalent sectors in the quantum theory. 

The integration over x is restricted to a finite interval Lq such that the Hamiltonian constraint appears simply as 


Ch 


N 


(fc2 2 fbSgn{p,) 

\/H 


+ 2bc^/\j^\ , 


( 3 ) 


in terms of rescaled variables given as 


b = b 

Pb = LoPb 


C = LqC, 
Pc = Pc- 


( 4 ) 
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Further, fixing the gauge N to j y^|pc| sgn{pc)/b, the Hamiltonian constraint can be further simplified to 

Ch = “+ 2cPc] ■ (5) 

With this simplified Hamiltonian and the Poisson bracket structure of the original phase space variables, one can obtain 
the equations of motion for these new rescaled variables whose solutions for exponentiated time variable (f = e^) 
appear as 

b{t) = — 1, 

c(i) = 

( 6 ) 


with m = to as an integration constant. 


The above solution gives the following line element corresponding to qab 


ds^ 



[vlf (^ 
Ll V t 


dx^ + t^dfl^. 


( 7 ) 


If one chooses the dimension of the integration cell Lq = this line element gives the interior of the Schwarzschild 
spacetime, and also the exterior when one moves pass t = 2m surface beyond which t appears as a spatial variable. 


B. Moving to the quantum corrected version 


For taking account of the quantum nature one then polymerizes the connections, i.e. uses the holonomies, 
of the connection A (parameterized by 5) instead of the connection itself and obtain one parameter family of the 
Hamiltonian constraint. The parameter 5 captures the notion of the length of the path along which the connection 
is integrated to obtain the corresponding holonomy. In this way the Hamiltonian constraint is also obtained up to 
the polymer parameter, in fact up to a parametric function A((5). All the elements of this family of parameterized 
constraints are required to reduce to the classical expression of the Hamiltonian constraint in the limit 5 —>■ 0. In 
terms of the canonical variable one can decide to make use of variable set sn^j of {b,c}. Thus, 

the Hamiltonian constraint in the gauge N = ( 7 -\/|pc|sgn(pc)^)/(sin A((5)(56) becomes 


Ca(s) = 


1 r sin(5c f smA{S)Sb \ t 

27 GArr 5 V ^ smA(S)Sb) 


( 8 ) 


Again, using the Hamiltonian constraint and the equations of motion one obtains 


where, 


, , 2 / ')5mp\ 

cyt) = - arctan =F 


PcW = ±^L 


2 t 2 


^ lompl y 


+ f 


cosA{S)5b = p{S) 


1 - 


f 2m 

V ‘ 


A{S)p{S) 


V{S) 


1 + 


f 2m 

V ‘ 


A{S)p{S) 


V{S) 


Pb{t) = - 


2 sin (5c sin A{6)6b Pc 
sin^ A((5)56 + 7^(5^ ’ 


p{S) = a/TTt^, 


ns) = 


+ ^ 2^2 _ 1 
y/l + + 1 


(9) 


( 10 ) 


Given these variables one can again construct the spacetime metric which will be bearing some implication of intro¬ 
ducing the polymerization in terms of dependence over b. Furthermore, as we can see that for d —>■ 0, everything 
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reduces to the previously discussed classical set-up, thereby assuring of the Schwarzschild limit when the polymer 
parameter S is insignificant. Therefore, in the case where the quantum correction is expected to be sub-dominant one 
is hopeful of obtaining the Schwarzschild metric back. Moreover, one also requires an asymptotically flat limit of this 
polymer corrected spacetime. One can verify that there exists only one such parametric function A((5) which ensures 
the asymptotic limit [9]. This choice for A((5) happens to be 


m 


1 

y/l + ' 


This modified spacetime has some salient features such as existence of two horizons, finiteness of the Kretschmann 
scalar with an asymptotic Schwarzschild behavior. For demanding mass independence of the location of Kretschmann 
scalar peak, one introduces a parameter oq of the full theory which is the minimum area gap in LQG. For a detailed 
discussions please refer to [9]. 

This approach of polymerizing the connection variable can in principle be implemented in two ways, namely where 
one polymerizes one of the component or both (b, c). The former is named as semi-polymeric treatment while the the 
second one is the full polymeric treatment. 

We consider the full polymeric loop black hole in this paper. This singularity free semiclassical black hole metric is 
given by 

ds^ = —f{r)dt^ + g{r)dr^ + h{r)r‘^did^, ( 11 ) 

where, 

(r-r+)(r-r_)(r-hn)2 
/W = - „4 I „2 -’ 


5(0 


(r — r+)(r — r_)r'^ 1 ^ 

{r + r*)2(r4 -|- ag) j 


where ag is the lowest area gap quanta in LQG, r_|_ = 2m, r_ = 2m'P^ and r* = ^r+r_ = 2mV, with the quantity 

^ V?TT-i 
v^^+l’ 

being the polymeric function which captures the correction from LQG in terms of e = jS. Qualitatively, this spacetime 
has similar horizon structure as that of a Reissner-Nordstrom black holes having two horizons at r_ and r+. Another 
important aspect of this black hole solution is that it is self-dual in sense of T-duality. One can verify that under the 
transformation r —>■ oo/r, the metric remains invariant, with suitable parameterization of other variables (details in 
[13]), hence marking itself as its dual description under T-duality . 

As we have argued before, the correction to observables because of e is of our interest and the correction due to 
oq is safely neglected; ag being a dimensionful quantity and hence appearing in observables with high energy scale 
suppression. So we will henceforth put ag to zero. Some interesting properties and phenomenology of this Loop 
Quantum Black Hole (LQBH) is explored in [13]. 


III. WEAK FIELD LENSING 


First let us consider lensing in the small deflection limit. If we expand the metric in the weak field limit, we get 

= + + ( 12 ) 


and 


g{r) = 1 -h ^(1 -H rf + ^(1 + 2V + 3ip2 -k 2V^), 


(13) 
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from which we can read off the ADM mass of the system M as 

M = 


(14) 


Rewriting the weak field expansion in terms of Arnowitt-Deser-Misner (ADM) mass and using the formalism of 
Keeton and Peeters [14] to identify the parametrized post-Newtonian (PPN) parameters, we arrive at the expansion 
form of the deflection angle in terms of the impact parameter b 


a{b) = Ai 




where the parameters in the expansion are given as 

Ai = 4 - 8 P and A 2 = 7r(f - ISV), 


(15) 


(16) 


keeping terms only up to first order in V since the polymeric function, 


_ \/TTe^-l 

V = 

y/TT^+l 


0{e^) < 1 . 


Precision solar system observations, more specifically solar gravitational deflection of Radio Waves using Geodetic 
Very-Long-Baseline Interferometry Data [15] put very strong constraints on Ai (the constraints on A 2 etc is very 
loose) as Ai = 3.99966 ±0.00090. Thus Ai is constrained to be extremely close to 4. This constraint can be expressed 
in terms of e (which we recall is the product of <5 and Immirzi parameter 7 ) and since the Immirzi parameter is 
estimated to be 7 ^ 0.25 by demanding that area law holds in loop quantum gravity (for example see [16, 17]) one 
gets 


(5<0.1. (17) 

It is surprising that solar system experiments have anything useful at all to say about a quantum gravity motivated 
parameter. This is because of a certain peculiarity of the model under consideration alluded to before (viz, the 
dimensionlessness of a quantum gravity motivated correction) which we discuss further later in this paper. We can 
also estimate the correction to various weak field observables (say the far field Einstein Ring (ER) location). Let 6^^ 
denote the angular location of ER for Schwarzschild metric. Then 






4Md;s 

doid 

OS 


where dis is distance between lens and source, dos is distance between observer and source and doi is distance between 
lens and observer. Using zeroth order solution to the lens equation following [14] we get the angular ER location for 
loop correction as 


Oe = OEoy ^ - &Eo{^ -'P), (18) 

up to first order in V. Thus ER shrinks slightly as compared to the Schwarzschild solution. 

There is a caveat however. We have assumed that the metric which was derived as a solution for black hole 
spacetime is a valid description for the vacuum spacetime geometry far away from a compact (horizonless) object in 
the theory. This is a non-trivial assumption but a plausible one. If these black holes can be formed via spherically 
symmetric gravitational collapse in the theory, starting from a compact object, the spacetime outside the object should 
not change in the process as the information about the collapse can not be communicated to asymptotia because of 
impossibility of monopolar gravitational wave generation. And so if the consequent black hole spacetime is unique 
the outside vacuum of a compact object must be described by the same spacetime asymptotically as that describing 
the black hole asymptotia. So with this working proposition, the equation Eq.(17) gives us a constraint on polymeric 
parameter S from solar system tests. 


IV. STRONG FIELD LENSING 


Now we turn to lensing in strong field lensing and explore the properties of relativistic images. 
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A. The formalism: strong lensing coefficients and observables 

It has been known since a long time that deflection angle of photons can become arbitrarily large (in principle 
infinite), if the photons are allowed to explore regions close enough to the black hole [18, 19]. So photons can loop 
around the black hole multiple times on the way from source to observer leading to the phenomenon of relativistic 
lensing and formation of relativistic images [20, 21]. The observation that the deflection angle diverges logarithmically 
as the distance of closest approach approaches the radius of the photon sphere has been used to define and compute 
the so called strong lensing coefficients and strong lensing observables by Bozza [22]. We use this formalism for our 
purposes. Since higher order images are clumped together, Bozza [22] considered a situation where one is able to 
resolve only the first relativistic image from the rest and compare the brightness of the first image with the cumulative 
brightness of all other relativistic images. The relavant strong lensing observables are then s, the separation between 
the outermost image and rest of the images; r, the luminosity ratio of outermost image and rest of the images (clumped 
together); and 0oo, the angular location of black hole shadow. Details of the calculation of strong lensing coefficients 
and observables is given in Appendix A. 


B. Strong lensing coefficients and observables for LQG Black Hole 

As we have seen earlier in Eq.(ll), in the Schwarzschild units the spacetime has double horizon at r_|_ = 1 and 
r_ = ■ This spacetime contains two photon spheres. The outer one is the unstable photon sphere while the inner 

horizon is the stable one. The radius of outer photon sphere (which is the one relevant for lensing) decreases with the 
increasing e as shown in Fig. 1. 



FIG. 1: The radius of photon sphere (scaled by Schwarzschild radius) as a function of e characterizing quantum gravity 
corrections 

As the radius of photon sphere decreases with the loop correction so does the corresponding impact parameter. 
Computation of the strong lensing coefficients shows that that ci and C 2 increase monotonically while Ups decreases 
monotonically as the quantum correction grows stronger (Fig. 2). 

At linear order in V, corrections to strong lensing coefficients can be calculated as 


Cl 


C2 


'^ps 


a.R 


4P f 2> 4P\ 

(1 + In (2 + T j + 

3^3 76r\ 


21 n( 6 ( 2 - ^)) + 


sVsr 

9(3 + V3) 


ln( 6 ( 2 - v^)) - 


8r[9 + V3 + ln C%\l'^ )] 

9(3 + 73) 


(19) 

( 20 ) 

( 21 ) 

( 22 ) 


As expected P —>■ 0 reproduces the Schwarzschild values. 

We now use these coefficients to compute the observational quantities r, 0^0 and s and plot their behavior with the 
increasing strength of the polymeric parameter in Fig. 3. The strong lensing observables have been computed for the 
galactic super-massive black hole scenario {M ~ 4 x lO®M 0 ,doj — 8.5kpc and with dos = 2doi ) as a concrete case. 
So, we see that with increasing polymeric parameter, the angular location of shadow region decreases slightly, the 
separation between the first and the last images increases, while the luminosity ratio decreases. The fractional change 
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FIG. 2: The strong lensing coefficients as a function of e characterizing quantum gravity corrections 




FIG. 3: The strong lensing observables as a function of e characterizing quantum gravity correction for the galactic super-massive 
black hole scenario, r is in magnitude and Ooo and s are in microarcseconds. 


is surely insignificant for practical purposes as for as one can foresee the limits of current and proposed observational 
technologies. However situation for time delay observations is more promising, as we shall discuss in next section. 
The time delay between successive relativistic Einstein rings is given by [23] (in Schwarzschild units, i.e., scaled by 
2M) 

. (jJr,:, I Ur,:, _E2_ / =JL\ , . 

Atat^at+i = 27rMp/i-b 2^.e ■=! (1-e‘^ij, (23) 

\ps y Cps \ / 

where 

2x - 1 

4x x^' 


c = 


(24) 





















To zeroth order time delay between successive relativistic images is simply 2'KUps x 2M and so as the critical impact 
parameter Ups decreases so does Ari 2 as shown in Fig. 4. 


Ar|_; 



FIG. 4: The time delay (in dimensionless units) between first and second relativistic Einstein rings a function of e characterizing 
quantum gravity correction 


V. OBSERVATIONAL PROSPECTS 

Since the factional change in angular location of photon sphere, which can be read off from Eq.(2I) is ^ the 
resolution needed for observing such a tiny change in relativistic images for i5 ^ 0.1 will be tens of nanoarcseconds as 
the typical angular location of the relativistic images in our galaxy is tens of microarcseconds 24/x"). This looks 
more optimistic if we consider billion solar mass black holes in other galactic centers as the radius of photon sphere 
will accordingly scale up. This will of course be compensated by distance factor to other galaxies when we consider 
the angular resolutions needed. However the time delay between relativistic images will be immune to the distance 
factor since time delays between successive images is roughly the is roughly the time taken by light to travel one 
closed orbit on the photon sphere. In fact it is well approximated by 2'iTUps [23] . So we can write 

AT = AT(1 - -^) (25) 

using Ups from Eq.(21) where AT is time delay between relativistic images for Schwarzschild black hole and AT is 
analogous quantity for loop black hole. The required time delay resolution to probe polymeric corrections for certain 
billion solar mass black holes can be of the order of seconds (assuming (5 ^ 0.1) as the time delay for them would be 
in hours [23]. 


VI. CONTRAST WITH OTHER QUANTUM GRAVITY SCENARIOS 

A notable feature of the quantum gravity motivated metric we study is that one of the parameters that characterize 
the deviation from Schwarzschild geometry is dimensionless. This is peculiar and is in contrast to other quantum 
gravity motivated scenarios like non-commutative black hole space time [24] where the non-commutativity parameter 
has dimension of or like the Callan-Myers-Perry black hole (CMP) [25] in String theory where the string tension 
has dimension of as well. In those scenarios non-commutative parameter and string tension set up a scale and all 
corrections are mass suppressed. In other words, if the quantum gravity motivated dimensionful parameter is Iqg then 

the corrections to observables generically takes the form where m is some positive power and M is the mass of 

the black hole. This mass suppression effect is absent in the space time under study and makes it more amenable to 
observational tests compared to other generic quantum gravity corrections such as the ones discussed before. Ignoring 
this mass suppression, however, strong lensing in non-commutative black hole space times [6] shows similar trend as 
the space time under consideration, and it might be interesting to investigate the effect of string tension corrections 
in CMP solution for strong lensing which is beyond the scope of current study. 


VH. EINSTEIN RING SYSTEM AS A PROBE: CONSISTENCY RELATIONS 

In a given model for space time describing a black hole one can look for useful relations between strong field and 
weak field observations and one can construct certain consistency relations. For example one can relate the far field 
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ER with the relativistic ERs (which are clumped near the angular location of photon sphere [22]). This was used in 
[26] to distinguish black hole and wormhole systems. This idea can be used to obtain consistency relation for strong 
and weak held lensing in observables for our model as below. 

Let Oq and be the far held and relativistic ERs for Schwarzschild spacetime and Oq and 0oo be corresponding 
quantities for the loop black hole. Then one can show that [26] 


Since we have 


eo = 9o{i-r), 


(27) 


from Eq.(18) and 


76P 

0oo=0oo(l-^), 


using the expression for Ups from Eq.(21) the consistency relation becomes 


- 373 ^ 0 ,- 2 / 58r\ 

V-—) 


(28) 


(29) 


for loop black hole. A deviation from this relation will invalidate the model. Since V is expected to be an universal 
constant, this gives a line in {Oq, 0 oo) plane where all observations must lie to be consistent with the model. 


VIII. DISCUSSION AND CONCLUSION 

The gravitational lensing signatures in both strong and weak dehection regimes were studied for some black hole 
solutions in LQG recently obtained by Modesto [9]. We highlighted a curious feature of the solution, viz., the 
appearance of a dimensionless parameter (related to Barberolmmirzi parameter and polymeric parameter) in the 
solution under study. This is in contrast to generic quantum gravity corrections such as the noncommutativity 
parameter or the string tension or the minimal area quanta etc. which are always dimensionful. This helps the 
solution to evade the mass suppression effects, i.e., the observables are corrected not by powers of ratio of quantum 
gravity scale to mass of black hole but by powers of the dimensionless parameter itself. By appealing to some 
general principles, we used the far away regions for this spacetime to model the spacetime away from any spherically 
symmetric matter distribution in the theory far away from central source. In particular then, this can be used 
to calculate corrections to observables in solar system observations. When used in conjunction with theoretically 
motivated values of Barbero-Immirzi parameter one can use this to constrain polymeric parameter from precision 
solar system observations. We find that the polymeric parameter S has to be less than ~ 0.1. We also computed the 
strong lensing observables and found that with increasing polymeric parameter, the angular location of shadow region 
decreases slightly, the separation between the first and the last images increases while the luminosity ratio decreases. 
Strong lensing coefficients are interesting even in themselves because of their relationship with quasinormal modes [7] 
and absorption cross sections from black holes [8]. Quasinormal modes have, for example, interesting link to semi- 
classical aspects of black hole physics [27] and for quantum gravity motivated solutions like the one studied in this 
paper. Therefore strong lensing coefficients might provide interesting theoretical motivations for further exploration. 

From observational vantage point, the magnitude of these effects is too small to be detected or ruled out in near 
future observations for galactic super-massive black holes. However, the required time delay resolution to probe 
polymeric corrections for certain billion solar mass black holes in distant galaxies can be of the order of seconds and 
seems a promising probe of these corrections once relativistic lensing observations become a reality. Even if we do 
not have the requisite angular resolution to distinguish between Schwarzschild and Loop Quantum Gravity motivated 
black holes, time delay observations might still be able to discriminate between these models. Important point to 
note is that this kind of quantum gravity correction has a reasonable chance of being probed in future astrophysical 
observations while noncommutativity parameter or string tension etc. have no realistic chance whatsoever because of 
the mass suppression effect. One caveat however is that any deviation from the vacuum solution in general relativity 
(for example black holes in Brans-Dicke theory etc.) would have qualitatively similar deviations from Schwarzschild 
strong and weak field lensing observables (although the exact strong-weak consistency relations would be quantitatively 
different) and one needs to be careful in interpreting the deviation from Schwarzschild spacetime. But the models 
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having the feature that the leading power of the parameters, in the expressions capturing the deviation (from their 
corresponding Schwarzschild values), is different in strong and weak field regimes (unlike Brans-Dicke where the 
correction shows up at linear order for both strong and weak field regimes) can be more easily distinguished from this 
LQG motivated model. 


Appendix A : BOZZA’S FORMALISM FOR COMPUTING STRONG TENSING OBSERVABLES 


The general spherically symmetric metric can be written in the form given by equation Eq.(ll) where f(r), g(r) 
and h[r) are arbitrary functions. We also restrict them to be such that the spacetime is asymptotically flat. 

For convenience we define two new functions in terms of metric coefficients 


w(r) 


Vf(.r)g{r) 

h{r) 


x{f) 


fir) 

Hr)' 


The photon sphere which describes a null geodesic at a fixed co-ordinate radius r is obtained by solving the equation 


rx'ir) = 2x 


(30) 


in terms of one of these new functions. The deflection angle which gives total angular deflection of a null geodesic 
coming from and going to infinity is given as 


Hro) = iHo) - TT, 


(31) 


where tq is the closest approach distance and 



Hr)r'^firo) 

Hro)ryir) 



(32) 


Let us define a dimensionless variable z = ^ and g = 1 — z. In terms of these new variables and the new functions 
described above, the expression Eq.(32) gets simplified to 


2 

We will be interested in the divergent structure of the integrand around r = tq which will give us a measure of the 
strong lensing coefficients, as we will see. Expanding the integrand in leading orders of g we separate out the integral 
into a divergent and a regular part as 


H ujjro/il-g)) 

Jo Vxiro) - xHo/H - ^))(1 - V)^ 


/(ro) = 2 


/ 


uj{ro/z) 


0 y/xHo) - xiro/z)z' 


zdz = 


with 


and 


/z3(ro) =2 [ 

Jo 


iHo) = IdHo) + IrHo), 

w(ro) 


0 \/7i(?'o)(l - z)+l 2 {rHi^ - z)' 


zdz, 


IrHo) = 2 [ 

Jo 


c 

uj{ro/z) 

w(ro) 

lo 

VxHo) - xirolz)z‘^ 

\/xiiro)H - ^) +72(»'o)(l - z)'^_ 


dz. 


In the above expressions 


(34) 


(35) 


(36) 


7i(ro) = 2x(ro) - rox'(ro), 


and 


l 2 {ro) = -xHo) + rox'iro) - {l/2)rlx''Ho). 
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The integral loiro) can be calculated exactly to give 

2w(ro) 

iDiro) = - , , In ■ 

V72(?’o) 


7i(?’o) 


y/li{ro) + l 2 (ro) + \/l 2 (ro) 


(37) 


which diverges for the tq = rps, with rpg being the radius of photon sphere. For rp close to Vps we can again expand 
loifo), we can write 


loiro) = -- 


V8uj{r 


"ps 


In I-1 I — In 2 


+ O{ro - rps), 


^2x{rps)-rl^X"{'rps) 

using 2x{rps) — rx'{rps) = 0. We also see that lR{ro) which converges for rp = rps, can be written as 

Iniro) = Inirps) + O{ro - rps). 

The deflection angle Eq.(31), which diverges for rp = rps, is thus given in the strong deflection limit by 


(r-p) = -oi In ( — - 1 ) + 02 + 0(rp - rps). 


' ps 


where 


(38) 


(39) 


Oi = 


V8uj{rph) 


l‘2.x(rps)-rlsX"(.rps) 
02 = -TT + an + qr, 
an = ui In 2, 


and 


aR = 


Inirps) = 2 /" 
Jo 



uj{rps/z) 

V^w(rps) 

lo 

Vxirps) - x{rpslz)z^ 

(1 - z)^2x{rps) - rl,x''{rps)_ 


dz. 


The deflection angle can be obtained as a function of the impact parameter u(rp) = 


\/x(7o) 


with 


where ci = ai/2 and 


a{u) = —Cl In (-1 ) + C 2 + 0{u — Ups), 

. ^ps 


ai , ^Xirps) - rl,x {rps) 

C2 = — In-^-l-a2- 

2 4x(rps) 


(40) 


(41) 


(42) 


where Ups = and doi gives the distance between the observer and the lens. In the above equation Eq.(42) 

\/x(rps) 

the coefficients Ci and C 2 are known as the strong leasing coefficients. They can be related to the leasing observables 
through the Virbhadra-Ellis lens equation [20] 


Further, we can express the deflection angle as a function of image position 6 using u = dpiO as 

a{9) = —Cl In (— -1^ + C 2 + 0{u — Ups), 


tan fd = tan 9 — C 3 [tan(Qf — 0) + tan 9] 


(43) 
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where C 3 = dis/dos- In the above expression jd gives the source position, while dis and dos are the lens to source and 
observer to source distances, respectively. For an aligned configuration we consider /3 and 9 to be small. Whereas, 
for strong lensing scenario the deflection angle a is typically of the order 2 n 7 r. Thus, n = 0 corresponds to the weak 
held limit which gives two weak field limit images while for n > 0 we get two infinite set of relativistic images. 

We substitute 


a = 2n7r + Aa^, 

with n G N and 0 < Aan -C 1, for one set of relativistic images (belonging to the same side of the source), the lens 
equation gets simplified to 


(3 = 9- c^Aan- 

The other set of images (on the other side of source) will then be given by 


(44) 


a = —2n7r — Aa„. 


Using Eq.(42) to the leading order we can invert a in favor of 9 as 


e{a) = ^ 1 + 

doi 


(45) 


We now make a first order Taylor expansion around a = 2 n 7 r, to express the angular position of the n-th image ss 

- CnAa„, (46) 

with 


/nO 


1 _|_ g(c2-2n7r)/ci 


and 


. _ Ups „(e2-2ra7r)/ci 

Sn — , c 

Cl djoi 


Solving the Lens equation upto first order gives. 


9n = el + ^{p-9<i). 

C3 


(47) 


(48) 


(49) 


For small ci and C 2 Eq.(48) suggests that is vanishingly small leading to the fact that all relativistic images lie 
very close to 0°. Similarly the other set of relativistic images have angular positions 


9n — —0° + — (/3 + 0))). 

C 3 

In the case of perfect alignment (/3 = 0), we obtain an infinite sequence of concentric rings, with angular radius 


(50) 


05 = ( 1 - ^ ] 00 
C 3 


(51) 


known as the Einstein rings (ERs). The magnification of the images is given by 


_ 1 

^”"/3 C3 


As we see from eqn Eq.(48) that (n gets exponentially suppressed for increasing n , the first image is the brightest and 
the brightness gets progressively decreased for increasing n. In accordance with Bozza [22] we consider a situation 
where one is able to resolve only the first relativistic image from the rest and compare the brightness of the first image 
with the cumulative brightness of all other relativistic images. Eor this purpose, we define lensing observables 


s — 91 9qqj 


(52) 
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the separation between the outermost image and rest of the images; and 


n—2 

the luminosity ratio of outermost image and rest of the images (damped together); 
as well as 


(53) 


^OO - ^ps/do/5 


(54) 


the angular location of black hole shadow. 

For high alignment approximation the above equations get simplified to 

s = (55) 

and 

r = - 1, (56) 

This gives a link between the calculated quantities ci,C 2 ,Ups and observational quantities r,s,0oo- 
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